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ABSTRACT 

The stellar winds flowing out of the components of WR + OB binaries 
can collide and the shock waves are formed. Stellar wind collision, particle 
acceleration by the shocks and generation of X-ray, 7-ray, radio and IR emission 
in WR + OB binaries are discussed. 
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1. INTRODUCTION 

One of the main properties of Wolf-Rayet (WR) stars is a very intense 

■ outflow of gas. The mass-loss rate for WR stars, M WR , and the terminal velocity 

of the matter outflow, V , far from the star amount to ~ (0.8 — 8) x 10 -5 M© 

yr" 1 and ~ (1 - 5)xl0 3 km s" 1 , respectively (Willis 1982; Abbott et al. 1986; 
t^. ' Torres, Conti and Massey 1986). No less than 40% of WR stars belong to 

binary systems. Young massive O and B stars are the secondary components 

^ [ of such systems. OB stars also have an intense stellar wind: M OB ~ 10 -6 M© 

yr -1 , V^ B ~ (1 - 3) x 10 3 km s -1 (Garmany and Conti 1984; Leitherer 1988). 

If the intensities of the stellar winds of WR and OB stars are more or less 
comparable or if the distance D between the components of the binary is large 



ON 



enough (see below), the winds flowing out of WR and OB stars can collide and 



Oh 



the shock waves are formed. In the shock the gas is heated to temperature ~ 10 



7 



K and generates X-ray emission (Prilutskii and Usov 1975, 1976; Cherepashchuk 
+3 . 1976; Cook, Fabian and Pringle 1979; Bianchi 1982; Luo, McCray and Mac Low 

1990; Usov 1990; Stevens, Blondin and Pollock 1992; Usov 1992; Myasnikov and 
Zhekov 1993) 

Until now X-ray emission from a few tens WR stars has been detected 
(Seward et al. 1979; Moffat et al. 1982; Caillaut et al. 1985; Pollock 1987). 
From the observational data it follows that the WR binaries are significantly 
brighter in X-rays than single stars (Pollock 1987), and the WR stars detected 
in X-rays have a much higher than usual incidence of binaries (Abbott and Conti 
1987). These two facts can be naturally explained if the X-ray emission of WR 
binaries is enhanced essentially by the radiation of the gas heated in the shock. 
Besides, the analysis of the energy spectra of the X-ray emission and the time 
variability indicates that the X-rays are generated in the outer regions of the 
stellar wind near the OB star (Moffat et al. 1982; Pollock 1987; Williams et al. 
1990), as expected in the model of stellar wind collision (see Section 2). 

Stellar wind collision may be responsible not only for the X-ray emission of 
WR + OB binaries and for their radio, IR and 7-ray emision as well (Williams, 
van der Hurcht and The 1987; Williams et al. 1990, 1992, 1994; Usov 1991; 
Eichler and Usov 1993). Below, the gas flow, stellar wind collision, particle 
acceleration, and generation of X-ray, 7-ray, radio and IR emission in WR + 
OB binaries are discussed. 



2. STELLAR WIND COLLISION IN WR + OB BINARY 



2.1 Stellar wind parameters 

For a typical WR + OB binary (M WR ~ 2 x lO- 5 M yr" 1 , M OB ~ 
lO~ 6 M0yr _1 , ~ V^ B ~ 2 x 10 3 kms _1 , the distance between the compo- 
nents of the binary D ~ 10 13 cm, the stellar wind gas temperature T~ 3xl0 4 
K) the parameters of the gas ahead of the shock fronts are the following: The gas 
density p~ 10 _14 g cm -3 , the sound speed V s ~ 10 km s _1 , the free path length 
of charged particles / ~ 10 2 cm, the Mach number £ = V fV 8 ~ 10 2 3> 1, and 
the Reynolds number 3? = (r OB /l) m ~ 10 12 . The free path length of particles 
and the Reynolds number behind the shoch fronts are ~ 10 9 cm and ~ 10 3 , 
respectively. The above values of 3? suggest that the gas flow is non-viscous and 
non-heat conductive. 

2.2. Geometry of the region of stellar wind collision 

The winds from the binary components are highly supersonic and flow 
nearly radially out to the shocks. In the shock the gas is heated to the temper- 
ature 

T(V ± ) = 2 7 ^4^m p Vl^3xl0 7 ( Vl - X K, (1) 
V ±J (7 + 1) 2 k p ± VlO 3 kms" 1 ; V ' 

where 7 is the ratio of heat capacities at constant pressure and at constant 
volume, n=A/(l + Z) is the mean molecular weight, A is the atomic weight of 
ion, Z is its electrical charge, k is the Boltzmann constant, m p is the proton mass, 
and V± is the component of the stellar-wind velocity perpendicular to the shock 
front. The value 7 is equal to 5/3 for a rarefied totally-ionized plasma. If helium 
predominates in the gas of the WR stellar wind, we have A = 4; Z = 2; \x = 4/3. 

Behind the shock the hot gas outflows from the region of stellar wind 
collision nearly along the contact surface (see in detail Usov 1992). 

In the case of the collision of two spherical winds with the terminal veloci- 
ties the distances r WR and r OB from WR or OB stars, respectively, to the region 
where these winds meet is 

"A r os = T ± Z^ D ( 2 ) 



(see Fig. 1), where 



= 1V± OB V OB 
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Since M OB ~ (0.01 — 0A)M WR and V^ B ~ the dimensionless param- 
eter r\ is small (rj <^1). Hence, the region of stellar wind collision is much nearer 
to the OB star than to the WR star (r WR 3> r OB ). 

The form of the contact surface C near the OB star (r ~ f OB ) is given by 
(see, e.g., Usov 1992) 

\Ux)\^r OB ^- x (3) 

for x ^ 7r /2 ( here x is the angle between the radius vector R c (x) horn the 
center of the OB star to the point at the contact surface and the line connecting 
the components of the binary). 



At intermediate distances (r OB Cr< (P/2)V^ R ) from the OB star the 
contact surface C approaches the conic surface C with angle 

9 ~ 2.1 (l - ^-^) r] 1 / 3 for 10" 4 < 77 < 1 . (4) 



The shock fronts £1 and £2 are near the conic surfaces Si and £2 at 
r OB < r < ( p / 2 ) V wr as welL In the case of wide WR+OB binaries, like WR 
140, when the fraction of the thermal energy of the hot gas lost by emission in 
the shock layers is small the angle A9 between these conic surfaces is ~ 9. In 
the other case the value of A6 is smaller than 9 and asymptotically approaches 
zero if this fraction goes to unity. 

The contact surface C and two shock fronts Si and £2 are twisted at a 
distance more than ~ [P/2)V^ due to the orbital rotation and have spiral form 
(see Fig. 1). Since the orbital velocities in massive binaries are typically much 
smaller than the terminal velocities we have (P/2)V°° ^> D. 

Since both r WR ^> r OB and the radius of OB star R OB is , as a rule, 
essentially larger than the radius of the WR star, in close binaries the gas velocity 
of OB star wind does not reach the terminal value before the shock S2. 

The velocity of the matter outflow V OB (r) at distance r from the center of 

the OB star varies from almost zero on the OB star surface to V for r > rter 
(Barlow 1982), where r ter is approximately equal to (3-5) R OB - 

If r OB > rter-, w e have the collision of two winds with terminal velocities 
and values r WR and r OB are determined by equation (2). If R OB < r OB < r ter , 
the gas flowing out of the OB star does not reach the terminal velocity by 
the time it enters the shock. In turn, the stellar-wind gas of the WR star is 
decelerated by the radiation pressure of the OB star while approaching the OB 
star. In this case, the gas velocities ahead of and near the shocks drop, and the 
gas temperature behind the shocks drops too (Usov 1990, 1992). If r OB < R OB , 
the OB wind is suppressed by the ram pressure of the WR wind, on the side 
facing the WR stars. 

2.3. Basic equations and boundary conditions 

The set of equations which describes the steady gas flow between the shock 
wave and the contact surface will be the continuity equation: 



the momentum equation: 
and the energy equation: 



div(pF)=0, (5) 
p(VV)V = -Vp , (6) 

P VV{H+\V\ 2 /2) = -Q , (7) 



where p is the density of the gas and Q is the energy loss per unit gas volume 
by radiation. 

Since the gas in the shock layer is practically totally ionized, its pressure 
p and its specific enthalpy H can be expressed as 



p= {N e + N l )kT= !— , 
rripp 



(8) 



7 p 5 fcT 
7 — 1 p 2 rripH 

here iVj = p/m p A is the ion density, N e = NiZ is the electron density. 

The ionized gas heated in the shock layer is emitting mainly due to free- 
free transitions of electrons in the Coulomb fields of ions at the gas temperature 
T > 3 x 10 6 K (Gaetz and Salpeter 1983 ). 

Gas parameters ahead of the shock front (index 1) and behind the shoch 
front (index 2) are related via the Rankin-Hugoniot relations 

P M n) =P*vi n \ Pi + Pi{vi n) f=p 2 + p 2 {V^f , 



V™=vi T \ H 1 + t{vi n) f = H 2 + t{V^f. (10) 

— * 

Indices n and r denote the normal and tangential components of the vector V. 
The condition V^ n) =0 is met on the contact surface. This condition and the 
Rankine-Hugoniot relations (10) are the full set of boundary conditions for the 
set of equation (5)-(9) needed to find the parameters of hot gas in the shock 
layers. 

2.4. The method of solution 

The set of equations (5)-(9) with the boundary conditions (10) was solved 
in (Galeev, Pilyugin and Usov 1989; Bairamov, Pilyugin and Usov 1990; Usov 
1992) analytically using the method of Chernyi (1961). In this method the ratio 
of the density of the gas, e, ahead of the shock to that behind it was considered 
as a small parameter, e < 1. If the fraction of the thermal energy of the hot gas 
lost by emission in the shock layer is small, e is equal to (7 — l)/(7 + 1) = 1/4. 
If the thermal energy losses are essential, the value of e decreases and goes to 
zero when practically all thermal energy of the hot gas is radiated. 

The problem of stellar wind collision was solved in the following way. 
First, the parameters of the hot gas in the external and internal shock layers 
were obtained as a function of the form of unknown contact surface. Then the 
form of the contact surface was found from the equality of the gas pressure on 
both sides of contact surface. 

In the Newtonian approximation, when we neglect the Busemann correc- 
tion, the equation for the contact surface is given by equation (3) . The difference 
between the analytical equation (3) and the numerical solution for the contact 
surface with Busemann correction is smaller than 10%. Therefore, to calculate 
both the parameters of the hot gas in the external and internal shock layers and 
X-ray emission from the shock layers equation (3) for the contact surface was 
used (Usov 1992). 

3. X-RAYS FROM THE REGION OF STELLAR WIND COLLISION 

3.1. Collision of two stellar winds with terminal velocities 

If a WR + OB binary is wide enough, r OB > r ter , the WR and OB winds 
collide with the terminal velocities, and the fraction of the thermal energy of 
the hot gas lost by emission in the shock layers is, as a rule, small. In this case, 



(9) 



the expected X-ray luminosity is (Usov 1992) 
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from the external shock layer between the shock front Si and the contact surface 
(see Fig. 1), and 



J int 



~ 1 3 x 10 35 ' MwR ' ' ^ OB ' ' WH 



lO^M©^" 1 i \ lO^Moyr" 1 / \10 3 /cms 
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^w^) (1555s) orgs " (12) 

from the internal shock layer between the shock front £2 and the contact surface. 

The total power of X-ray emission from the region of stellar wind collision 
is L x = L ext + Li nt . The X-ray spectra of the shock layers may be roughly 
approximated by the bremsstrahlung radiation of the uniform totally ionized 
plasma with the temperature ~ (0.8 — 0.9)T(V™ R ) for the external shock layer 
and ~ (0.8— 0.9)T(V^ B ) for the internal shock layer, where T is given by equation 
(I)- 

To calculate the parameters of hot gas in the shock layers and the X-ray 
emission of this gas the values of e and ry 1 / 2 were considered as small parameters. 
Since the zero and first terms of the series were taken into account only, the 
accuracy of our calculations is of the order of (e 2 + if) 1 / 2 , i.e. the accuracy is 
several tens of percent for 77 < 0.1 which is typical for WR+OB binaries. 

In our consideration of X-ray emission from the region of stellar wind 
collision we did not take into account the absorption of X-rays in the stellar 
winds. We can do it for wide WR + OB binaries, D > 10 13 cm, (Usov 1992) 
when the sightline to the region of stellar wind collision is passed far enough 
from the WR and OB stars. 

3.2. Close and very wide binaries 

In a close binary, r OB < r ter , both the WR stellar wind and the OB stellar 
wind have velocities near the shock wave which are less than their terminal 
velocities (see Section 2.2). To estimate the X-ray emission from the region of 

stellar wind collision it is necessary to substitute V and V" for V* and 

J WR OB WR 

V* B in equations (11) and (12) (here V£ and V* B are the gas velocities of the 
WR and OB stellar winds, respectively, ahead of and near the shocks). Let us 
discuss briefly the behaviour of the X-ray luminosity, L x , for a remote observer 
when the value of r OB varies from r ter to R OB . 

At first the X-ray luminosity L x may increase a few times in comparison 
with the sum of equations (11) and (12), L x oc (V^^) -1 . In turn, the gas 



temperature in the shock layers drops, T oc (V^. OB ) 2 , and the spectrum of the 

X-ray emission shifts to the soft region (Usov 1990). Then, if the value of r OB 
is a few times smaller than r ier , the main portion of this soft X-ray emission is 
absorbed in the WR and OB stellar winds and the X-ray luminosity is decreased 
(Luo, McCray and Low 1990; Usov 1990). At 1 < (r OB /R OB ) < 1.2-1.3 the X- 
ray emission is suppressed practically at all. Hence, it is expected that the X-ray 
luminosity of a WR + OB binary is as high as possible when r OB ~ (0.5 — l)r ter 

orD~ (0.5 - l)??- 1 / 2 ^. 

It is worth noting that some increase of the X-ray flux may be observed 

near the conic surface C from more or less close WR + OB binaries for which the 
X-ray absorption in the stellar winds is essential. The point is that the outflowing 
plasma in the shock layers is highly ionized, and the X-ray absorption along this 
surface may be small. 

Equations (11) and (12) are valid only if there is temperature equalization 
between ions and electrons in the shock layers. If there is no such an equalization, 
the power of the bremsstrahlung emission from the layer is ~ L(D/.D e(J ) _1 / 3 at 
D > D eq , where L is given by either equation (11) or equation (12) depending on 
the kind of the shock layers and D eq is the distance between the conponents of 
the binary at which temperature equalization in the shock layer is broken (Usov 
1992). The temperature of electrons which determines the X-ray spectrum is 
oc (D/D eq )~ 2 / 3 at D > D eq . 

Among WR + OB binaries, WR 140 is the most powerful X-ray source, 
L x ~ 4 x 10 34 ergs s _1 (Pollock 1987). It is natural because at the moment of 
the X-ray observation WR 140 was more or less near the state with r OB ~ n er . 
The observed value of L x is in agreement with the X-ray luminosity which is 
expected from the region of stellar wind collision in WR 140 (Usov 1992). 

Except for the X-ray emission generated in the region of stellar wind colli- 
sion, WR + OB binaries may have an X-ray emission which is inherent in single 
stars. X-ray emission of single massive stars may be generated due to heating 
of the outflowing gas to temperatures of about a few million degrees either by a 
large number of quasi-periodic strong shocks (Lucy and White 1980; Cassinelli 
and Swank 1983; Owocki and Rybicki 1985) or by the current sheets (Usov and 
Melrose 1992). From analysis of time variability of X-ray emission from WR 
+ OB binaries it is possible to separate X-ray emission of colliding winds from 
X-ray emission that is inherent in single stars. 

4. GAS COOLING AND DUST FORMATION 

In the shocks the gas of stellar winds is compressed but not more than 
~ 4 times. If a WR + OB binary is close enough and during the gas motion 
through the shock layer the main part of the thermal energy of the hot gas 
is radiated, the degree of the gas compression can increase significantly, i.e. a 
cold dense gas may be formed. The terminal velocity of the cold gas has to 
be essentially smaller than V WR and V QB . This cold gas may be heated and 
accelerated because of its interaction with the WR high-velocity wind. 

In the case of long-period binaries like WR 140 even near the periastron 
passage the fraction of the thermal energy of the hot gas lost by emission is 
small. At first sight in this case the hot gas outflowing in the shock layer has to 
expand quasi-adiabatically, and the essential cooling and compression of the gas 
in the region of the stellar wind collision is impossible. But it is not so. Indeed 



the main part of the hot gas in the shock layer expands quasi-adiabatically. 
However, near the contact surface there is the region where the energy losses 
may be important and the gas may cool and compress very strongly (Usov 1991). 
The strong cooling in this region is caused by the slow motion of the gas along 
the contact surface, such that the gas has ample time to cool. In this case the 
fraction of the WR wind which may be strongly cooled and compressed in the 
shock layers is (Usov 1991) 

a ~ r?( 2 (-?-) " ( V - 
1 yiO-SMoyr- 1 ) \10 13 cm J \1Q* cms- 1 




The characteristic cooling time of the hot plasma via the bremsstrahlung 
depends on the temperature and the density as r c ~ 1.5 x 10 11 A^ _1 T 1 / 2 s. Taking 
into account that the external pressure of the hot gas compresses the cooling 
gas (N oc T _1 ), we have r c oc T 3 / 2 , i. e. the process of cooling near the contact 
surface is accelerating (thermal instability). The gas near the contact surface 
may be cooled down to ~ 10 4 — 10 5 K, and the density of the cold gas may be 
as large as ~ 10 3 — 10 4 p*, where p* is the gas density before the shock (Usov 
1991). 

At the distance r<i — a few xlO 15 cm from the binary, the cold dense gas 
can form a dust shell and reradiate the UV radiation in the IR region (Williams 
et al. 1990, 1992, 1994; Usov 1991). 

The value of a is, as a rule, very small at r OB > rter because of high 
velosity of the stellar winds and increases very sharply at r OB < r ter - When 
r OB is comparable with R OB , the dynamical gas pressures of both winds ahead 
of the shocks are decreased, and the gas compression in the shock layers is 
decreased too. In addition, in close WR + OB binaries with r OB ~ R OB the 
X-ray absorption inside the shock layers is high, and it can prevent a strong gas 
cooling in the layers. Therefore, when the value of r OB is somewhere between 
~ r ter and ~ 0.3rt er (say, r OB ~ 0.5n er ), the most favourable conditions for a 
strong compression of the gas in the shock layers are realized (Usov 1991). 

The long-period WR binaries (WR 48a, WR 125, WR 137, WR 140, and 
possibly WR 19) with a large enough value of the orbital eccentricity are dis- 
tinguished among the rest of the WR binaries in that they experience the state 
with r OB ~ 0.5rter in the process of their orbital motion. In this state there is a 
strong outflow of cold gas, which leads to an IR outburst (Williams et al. 1992, 
1994). 



5. PARTICLE ACCELERATION AND RADIO AND 7-RAY EMISSION 

It is well known that plasma shocks in astrophysical setting can and do 
accelerate charged particles to high energies (for a review, see Blandford and 
Eichler 1987 and Ellison and Jones 1991 and references therein). Therefore the 
strong shocks formed by the colliding stellar winds in WR + OB binaries are 
very promising as particle accelerators. 

Necessary conditions for particle acceleration by a shock wave include the 
following (Eichler and Usov 1993): 

1. The shock must be collisionless in the absence of an external injection 
mechanism. This is typically satisfied if the strength of the magnetic field ahead 



of the shock is B ^> 10 -6 Gauss. Here and below all numerical estimates are 
given for a typical WR + OB binary (see Section 2.1). 

2. The Ohmic damping rate of Alfven waves must not exceed the maximum 
growth rate. 

3. For primary electron acceleration, the shock velocity must considerably 
exceed the phase velocity of whistler waves propagating normal to the shock. 

Last two conditions holds if B <C a few Gauss. 

All of these constraints on B are typically satisfied by colliding winds in 
WR + OB binaries. Hence the region of stellar wind collision in these binaries 
may be a strong source both high energy particles and nonthermal radiation 
generated by these particles. 

The maximum electron energy is limited by ion-neutral wave damping and 
by inverse Compton losses (Bell 1978). The maximum energy allowed by ion 
neutral damping should be well above 100 Gev and does not play an important 
role. 

The upper limit on the Lorentz factor of electrons due to inverse Compton 
losses can be shown, assuming a quasi-parallel shock geometry, to be given by 
(Eichler and Usov 1993) 

imax " Xa T L bol { c J ~ 6XLU ^ Gauss J ' (A4J 

with approximate equality obtaining in the absence of any other, stronger limits. 
Here L^i ~ 10 39 ergs s _1 is the bolometric luminosity of the OB star, a T is the 
cross-section of Thomson scattering, e is the electron charge and A ~ 3 is the 
ratio of mean free path to gyroradius. 

Given the constraints of B discussed above, the energy of electrons i? max = 
r max mc 2 is enough to generate synchrotron radio emission at the frequencies of 
which nonthermal radio emission from WR + OB binaries has been observed 
(Florkowski and Gottesman 1977; Becker and White 1985; Abbott et al. 1986; 
Felli and Massi 1991). 

A differential energy spectrum expected for relativistic electrons acceler- 
ated in the region of stellar wind collision is N(E) oc E~@ , where N(E)dE is 
the number of high energy electrons per unit volume with total energy E in the 
interval dE and (3 ~ 2. The spectrum of synchrotron radiation generated by 
relativistic electrons with this energy spectrum is S(y) oc u~ J , here S(u) is the 
flux of radiation at the frequency v and 7 = [fi — l)/2 ~ 0.5 is a spectral index. 

Typically, only a small fraction of the energy of accelerated electrons are 
converted to radio emission before the flow convects them out of the colliding 
wind region. 

From equation (3) we can get the characteristic size of nonthermal radio 
emisson 

l rad ~ 2R c (n/2) ~ nr OB . (15) 

Stellar winds outflowing from the WR and OB stars are opaque to free- 
free absorption. A typical radius of the radio photosphere at the frequency of 
a few GHZ for WR stars is R rad ~ a few times 10 14 cm (Wright and Barlow 

1975). The value of this radius for OB stars, R™ B , is approximately an order 
of magnitude smaller than R rad . Hence, in the case of WR and OB stars the 
radius of their radio photosphere is of the order of hundred or more stellar radii. 



Radio emission from the colliding wind region may be observed at the 
frequency v only if the time R™n/V™ R during which the gas outflowing from 
the WR star with the velocity reaches the radius of the radio photosphere 



R™d j s of the order of or smaller than the time which is necessary for the 



W R 



WR+OB binary to turn on the angle 29 + A9 ~ 39 (see Fig. 1), i.e. 

r>rad no n r>rad —2/3 

For a typical WR+OB binary (see Section 2.1) the value of P cr is of the 
order of a month when the frequency of observation is a few GHz. 

If P < P cr nonthermal radio emission from the region of stellar wind 
collision in the WR+OB binary cannot be observed because it is screened by 
the WR wind. 

Striking variations of radio emission from remarkable WR binary WR 140 
were observed (Florkowski and Gottesman 1977; Becker and White 1985; Abbott 
et al. 1986). This time variability and the other observational data on radio 
emission of WR 140, such as luminosity, spectrum, the size of nonthermal radio 
source, etc., may be explained by a model in which its nonthermal radio emission 
is generated at the site of stellar wind collision (Eichler and Usov 1993). 

Gamma-ray emission can be attributed to either inverse Compton emission 
or p — p collisions followed by pion decay. In the former case, the seed photons 
would presumably be the processed UV radiation from the OB stars. A high- 
energy tail of the inverse Compton spectrum of 7-rays may exist up to ~ 10 3 
MeV (Eichler and Usov 1993). 

Considerations of shock acceleration theory suggest that the WR 140 sys- 
tem may be a strong source of 7-rays (Eichler and Usov 1993). The 7-ray flux 
has to increase essentially near the periastron passage. Recently, the hard X- 
ray and soft 7-ray emission from WR 140 was observed by the OSSE detector 
aboard the CGRO (Hermsen et al. 1994). 
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